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Abstract - In this Letter, we present the first experimental study of bridge structures in three- 
dimensional dry granular packings. Our experiments, which use X-ray tomography, are in agree- 
ment with both present and earlier simulations. When bridges are small, they are predominantly 
'linear', and have an exponential size distribution. Larger, predominantly 'complex' bridges, are 
confirmed to follow a power-law size distribution. Simulation and experimental results for the 
distribution of sizes, end-to-end lengths, base extensions and orientations of bridges show good 
qualitative and quantitative agreement. This robustness underlines the fundamental role played 
by bridges as spatial heterogeneities in granular media. 
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The study of random granular packings remains an ac- 
tive research field [T]. For packings containing frictional 
grains, it is now well established that cooperative struc- 
tures such as bridges, are ubiquitous: these are defined as 
collective structures where neighboring grains rely on each 
other for mutual stability [2!. Apart from their intrinsic 
interest, bridges are intimately related to the phenomenon 
of compaction [3^ ; it has been shown that much of the com- 
paction to high densities in shaken packings occurs via the 
cooperative dynamics of bridge collapse [4] . Bridges can 
be further classified [5^ as linear or complex, depending on 
the topology of their backbones or contact networks: com- 
plex bridges have backbones with loops and/or branches, 
while linear bridges do not. Their dynamical implications 
are equally interesting: it has been suggested that grain 
motion close to the jamming transition is via the mo- 
tion of 'dynamic linear chains' [lllH], which are akin to lin- 
ear bridges Force networks in anisotropically sheared 
static packings [9l[T0] show fractal dimensions similar to 
those of linear bridges formed under gravity, thus rein- 
forcing connections that have already been made j3| [TT| rr2] 
between linear bridges and force chains. 

Recent experimental advances have allowed for the non- 
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invasive imaging of structure in dry granular media, using 
for example magnetic resonance imaging or positron emis- 
sion tomography [131 [TJ]. Although interesting attempts 
have been made to use these tools to characterise force 
chain distributions, they rely on indirect geometric mea- 
sures of, say force chain lengths [15], rather than direct 
measurements of force. Bridge structures have also re- 
cently been probed experimentally in colloidal packings 
p!6] and compared with the results of computer simula- 
tions of shaken (dry) granular media [5]; while the exper- 
imental results are of interest in and of themselves, it has 
to be remembered that colloids are governed by thermal 
energy, while temperature does not govern the dynamics 
of dry granular media. Instead, a perturbation such as 
shaking or tapping has to be applied [17, in order to gen- 
erate particulate motion in dry granular packings; such 
athermal perturbation could conceivably generate its own 
particularities in bridge structure, absent in the colloidal 
case. It is therefore important to characterise the statis- 
tics of bridges in the steady state of shaken dry granular 
packings in 3d and to compare them with the results of 
computer simulations modelling exactly the same physical 
situation; this is the main purposes of this Letter, where 
we use direct and non-invasive tomographic measurements 
to characterise bridge structures in dry granular packings. 

In the experiments, two packings (one monodisperse 
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(Duke Scientific, USA) and one polydisperse) of glass 
beads were used, with packing fractions of 0.623 and 0.597 
respectively. The glass beads had diameters of 200 ±15 fim 
and 300 ±50 fim respectively: there were ~ 18000 beads in 
the monodisperse packing and ~ 5000 beads in the poly- 
disperse packing. Both packings were tapped for more 
than 10, 000 cycles using a commercial shaker to ensure 
that a steady state had been reached. The tapping proto- 
col involved a single 30-Hz sine wave at a rate of 1 Hz with 
an effective tapping amplitude of 2.82 g. An X-ray micro- 
tomography machine (MicroXCT-200, Xradia Inc., USA) 
was used to make structural measurements, with 1200 pro- 
jection images taken on the samples. The tomography- 
reconstructed 3d images were analyzed by a marker-based 
watershed imaging segmentation technique p!8l[T9] using a 
custom-made Matlab program. A series of preprocessing 
procedures were adopted to remove artifacts and enhance 
the images, as shown in Fig. [TJ thus allowing for accu- 
rate identification of grain sizes and positions. The high 
spatial resolution facilitated the accurate determination 
of contacting neighbors, crucial for the identification of 
bridges. 
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Fig. 1: (color online) Image processing steps from the raw re- 
constructed image to the segmented image using a marker- 
based watershed algorithm. 
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Fig. 2: (color online) (a) Log-linear plot of experimental size 
distributions P{n)iin of linear bridges: the full and dashed 
lines represent exponential- law fits for mono- and polydisperse 
beads respectively. The corresponding exponent values are 
a = 0.45 ± 0.02 (monodisperse) and a = 0.49 ± 0.04 (poly- 
disperse). (b) Log- log plot of experimental size distributions 
P{n)comp for all bridges. The power-law fits to the largely com- 
plex bridges for n > 7 yield r = 2.1 ± 0.2 (monodisperse) and 
T = 1.9di0.4 (polydisperse). (c) Log-linear plot of size distribu- 
tions P(n)iin for linear bridges in packings at different densities 
(j) from computer simulations of shaken monodisperse spheres. 
The black line is an exponential- law fit as in panel (a) to the re- 
sults for (j) = 0.62, yielding a = 0.46 ±0.01. For the three lower 
densities, data collapse is obtained, giving a = 0.49 ±0.02. (d) 
Log- log plot of size distributions P{n)comp for all bridges in 
packings at different densities 0, from computer simulations 
of shaken monodisperse spheres. The black line represents a 
power-law fit for n > 7 as in panel (b) for (j) = 0.62, giving 
T = 2.6 ± 0.1. The exponent r appears to vary strongly with 
packing fraction. 



The identification of bridges in experiment and simu- 
lation followed the overall algorithm used in [5l|Tl][20]. 
Sphere coordinates are transformed into a list of at least 
three contacts for each particle: each particle in a contact 
list is next identified with a unique set of three other par- 
ticles that provide its supporting base, using in this case 
a criterion favoring the stabilizing triplet with the lowest 
mean-squared separation [21]. One then looks for sets of 
particles that appear in each others stabilizing triplets : 
these are clearly particles which are mutually stabilizing, 
and are thus identified with bridge particles. Finally, clus- 
ters of mutual stabilizations can be identified to reveal a 
unique set of bridges in a static close packing. 

Our experimental findings suggest that linear bridges 
predominate for sizes of up to n 10, which are char- 
acterised by a simple exponential distribution P(n) 

lin 

g-an (Fig. [2ja)). For larger sizes n, complex bridges pre- 
dominate, which are in turn characterised by a power-law 
distribution P{n)comp ^ (Fig. [2fb)). These results 
are robust to the presence of polydispersity: remarkably. 



even the associated exponents agree - within error bars - 
in the two cases, with a = 0.45 ± 0.02 and r = 2.1 ± 0.2 
(monodisperse) and a = 0.49 ± 0.04 and r = 1.9 ± 0.4 
(polydisperse). 

The above results are in broad qualitative agreement 
with earlier simulation results on simple and complex 
bridges [5 , which is in itself quite good. There are, how- 
ever, quantitative differences in the values of the expo- 
nents measured. We speculate that this could be due to 
the difference in packing fractions (j) considered in the ear- 
lier simulations and in the current experiment (0.56 and 
0.62 respectively); since packing fraction is a fundamental 
structural descriptor of granular media, it could reason- 
ably lead to quantitative, if not qualitative differences in 
the size distribution of granular bridges. The simulations 
described below investigate this issue, and confirm this 
dependence. 

Accordingly, we generated configurations corresponding 
to (j) = 0.58, 0.59, 0.60 and 0.62 using a well-established 
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hybrid Monte Carlo sphere shaking algorithm [22]. The 
simulations were performed on 1630 spheres in a rectan- 
gular cell with lateral periodic boundaries and a hard dis- 
ordered base, using 100 different random initial configura- 
tions per shaking amplitude. Care was taken to use stable 
configurations in the steady state, saving about 200 stable 
configurations (picked out every 500 cycles to avoid corre- 
lation effects) for bridge identification and analysis. Our 
findings are shown in Figs. [2] (c) and (d) where we find 
that there continues to be excellent qualitative agreement 
between experiment and simulation. Quantitatively, we 
observe an interesting trend: the agreement between ex- 
periment and simulation gets better as the corresponding 
packing fractions converge. In fact for (f) = 0.62, our simu- 
lations yield values of a = 0.46±0.01 and r = 2.6±0.1 (see 
black line in Figs. [2] (c) and (d)). The perfect agreement 
within error bars for the a values between experiments 
and simulations conducted at the same packing fraction 
suggests that this is the most important parameter con- 
trolling the behaviour of linear bridges; the very slight de- 
viation for r suggests that another descriptor, such as the 
coordination number, could play an important secondary 
role in the case of complex bridges. We notice also that 
while the results for the three lower packing fractions in 
Figs. [2] (c) and (d) can be collapsed onto the same straight 
line (yielding a common slope and hence a common expo- 
nent), those for (j) = 0.62 are significantly different, and 
need to be plotted separately. It is tempting to specu- 
late that the onset of ordering (in both experiment and 
simulation) is responsible for this quantitative change in 
exponent values ^. 

All these trends persist in measurements of the other 
bridge structure descriptors considered below; in the fol- 
lowing we focus on linear bridges, leaving the detailed ex- 
amination of complex bridges to future work. The first of 
our descriptors is the base extension, whose definition we 
review. If all possible connected triplets of base particles 
for a particular bridge are considered, the vector sum of 
their normals is defined to be the direction of the main 
axis of the bridge, typically inclined at some angle 6 (the 
orientation angle) to the z-axis. The base extension b is 
defined as the radius of gyration of the base-particles about 
the z-axis^ and is a measure of the spanning or jamming 
potential of a bridge [5]. 

In Fig. [3l we present experimental and simulation re- 
sults obtained for the base extension distributions of lin- 
ear bridges. Figs. [3] (a) and (c) correspond respectively to 
experimental and simulation plots of p{b\n) (normalised 
probability distributions of b conditional on bridge size n) 
vs. b for different bridge sizes n. They are qualitatively 
similar, showing sharply peaked distributions which flat- 
ten out with increasing bridge sizes, as seen in earlier sim- 
ulations [5 . Figs. [3] (b) and (d) correspond respectively 
to experimental and simulation plots of log p{b) (cumu- 
lative probability distribution of b) vs. the normalised 
variable b/ < 6 >, with < b > the mean extension of 
bridge bases. Both show the exponential tail in the dis- 
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Fig. 3: (color online) Probability distributions, for monodis- 
perse beads, of base extensions b for linear bridges, (a) Exper- 
imental plot of p(6|n) vs. 6 for different n. (b) Experimental 
plot of the log of the cumulative distribution, log {p{b)) vs. 
b/ < b >. (c) Plot of p{b\n) vs. b for different n, generated 
from computer simulations of shaken monodisperse spheres, 
(d) Plot of log {p{b)) vs. b/ < b > from computer simulations 
of monodisperse sphere packings at different densities. 



tribution function noted in the results of earlier simula- 
tions [5], suggesting that bridges with small base exten- 
sions are not favored. This result appears to be robust 
both with respect to polydispersity in the experimental 
results and packing fractions in the simulations. It also 
reinforces earlier suggestions [3l of deep connections be- 
tween force chains and linear bridges: the cumulative dis- 
tributions of force chains in anisotropically sheared gran- 
ular systems [T0l[23] as well as MD simulations of anal- 
ogous particle packings [231 [25] show very similar expo- 
nential tails. Recent simulations have directly confirmed 
the connection between force chains and linear bridges, 
by suggesting that forces are principally transmitted by 
particles in bridges [12]. 

Another important quantity related to a linear bridge 
of size n is obviously its 'span', i.e. its rms end-to-end 
length Rn. Our results for this are presented in Fig. Sj^a) 
(experiment) and Fig. |4]^b) (simulations). As expected, we 
find the scaling law Rn ~ in both experiment and 
simulation. We obtain exponent values u = 0.61 ±0.02 
from experiment and u = 0.60 ±0.01 from simulations. 
Again, the agreement between experiment and simulations 
is remarkable; additionally, we do not observe a strong 
dependence on packing fraction. Given that 0.59 for a 
3d self-avoiding random walk, this suggests that the linear 
bridges that we have examined here look - within error 
bars - exactly like self-avoiding walks in three dimensions. 

Finally, we examine via experiment and simulation, the 
normalised distribution for the mean angle made by a 
linear bridge with the z-axis, as well as that of its variance 
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Fig. 4: (color online) Plot of the (log of the) end-to-end length, 
Ini^n vs. ln(n — 1) for linear bridges in monodisperse sphere 
packing. The full lines represent power-law fits. For the exper- 
imental plot (a) V — 0.61 di 0.02, while (b) simulations for all 
the densities give v — 0.60 ± 0.01. 



[5] in Fig. [5l In Figs. [5] (a) and (c), we plot respectively 
experimental and numerical results for p(0|n), the ori- 
ent at ional distribution conditional on n. The cumulative 
distributions p(6), also plotted in these figures, closely re- 
semble each other as well as confirming earlier results 
Notable features are a peak around 20^ as well as the de- 
crease of with increasing bridge size, which [5^ suggests 
that larger linear bridges form domes. In Figs. [5] (b) and 
(d) we plot the variance of the mean orientational angle 
(6^), against size n. According to [5 , (B^) obeys 
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where is the equilibrium value of the variance of the 
link angle. This is plotted as the full line in the figures, 
where the symbols represent results for experiment and 
simulation respectively. These look similar to each other, 
and also to the results of earlier simulations [5l. 

We now discuss our results. The experiments reported 
in this Letter show satisfying qualitative agreement with 
the results of present and earlier [5] simulations, in the 
sense that the forms of distributions of quantities ranging 
from sizes to orientations and base extensions are robustly 
the same in every case. This is already quite good for such 
a complex field, and suggests that there is truth in the 
idea of bridges being good characterisers of spatial het- 
erogeneities in granular media. Quantitative differences, 
where these exist, between the current experiment and 
earlier simulations [5] have been successfully ascribed to 
the fact that the data were taken at rather different pack- 
ing fractions 0: the simulations reported in this paper have 
been conducted at a range of packing fractions and man- 
ifest this dependence, such that there is excellent agree- 
ment between simulation and experiment at matching 0. 
This reinforces the intuitive idea that packing fractions 
should influence the details of bridge structure (while leav- 
ing the qualitative features unchanged); higher packing 
fractions, for instance, probably constrain linear bridges 
to be fully three-dimensional rather than leaving them the 
choice to be planar. We might expect that for complex 
bridges, the factors at play (apart from the robust power- 
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Fig. 5: (color online) Orientational distributions of linear 
bridges in monodisperse systems from experiment and simu- 
lation, (a) Plot of the orientational distribution conditional 
on size n, p(B|n) vs. G, obtained experimentally, (b) Experi- 
mental plot of the variance of G, (G^^ as a function of n\ the 
full line shows the fit to Eq. [T] (c) Plot of the orientational 
distribution conditional on size n, p(G|n) vs. G, obtained from 
computer simulations of shaken spheres, (d) Computer simu- 
lation generated plot of the variance of G, (G^) as a function 
of n\ the full line shows the fit to Eq. [T] as in (b). 



law size distribution) may indeed be more complex: since 
these are branched structures, we might expect that co- 
ordination number, might have an important role to play 
in precise quantitative measurements of structure, as in- 
dicated by our measurements of the exponent r. We hope 
to carry out a detailed study of complex bridges in future 
work. 
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